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We provide an analytic solution to the general wavelength integro-differential equation describing 
CN| ■ the damping of tensor modes of gravitational waves due to free streaming neutrinos in the early 

universe. Our result is expressed as a series of spherical Bessel functions whose coefficients are 
U . functions of the reduced wave number Q. 
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INTRODUCTION 



Observations of the cosmic microwave background (CMB) have given increasing support to inflationary cosmo- 
logical models. Density perturbations during this inflationary period are believed to have given rise to the large 
[ scale structures of the universe [![. In addition to these scalar perturbations, a spectrum of gravitational waves 
is also produced Q (tensor modes) which could provide information about the early universe. In particular, 
i ' the contribution of these tensor modes (measured in terms of a tensor-to-scalar ratio Q) to the temperature 
anisotropy of the CMB could be used to check the predictions of inflationary models. Here we shall consider- 
ed " the effect of anisotropic inertia on the gravitational radiation and confirm that it is non- negligible. Although it 
will be assumed here that the anisotropic inertia is dominated by neutrinos and antineutrinos Q, Ref. Q has 
proposed a free streaming relativistic gas contribution and a method to constrain its fraction density through 
\ measurements of the CMB B-polarization spectra. 

r***- \ In a 2004 paper [|[ , Weinberg derives an integro-differential equation for the propagation of cosmological 

j | . gravitational waves. He writes a wave equation for the perturbation to the metric hij(x.,t) and defines x( u ) as 

fctf=M0)x(«), (!) 



where u is the conformal time multiplied by the wave number 

rt dt' 



a(t>) 



(2) 



The wave equation for the perturbation leads to an integro-differential equation which the function \ satisfies for 
general wavelengths. In the variable y = a(t)/aEQ, where cieq is the expansion parameter at matter-radiation 
equality, this is Eq. (32) of 



/ 2(l+j/) 1\ , Mn 2 , , 24/„(0) f y ,Mrf) A , 

0- + y)x (y) + I — - — + 2 ) ^ yy> + Q xw = — ~^ — J K (y>y) dyl d v ■ 



(3) 



with the initial conditions: 



X (0) = l, X'(0) = 0. (4) 

Here /„(0) = 0.40523 is the fraction of the energy density in neutrinos and Q = y/2k/kEQ- The kernel K will 
be discussed below, and y is related to u in the following manner 



2Q Vl + y-l . (5) 
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We will provide an analytic solution of Eqs. ([3]) and (U]) that is valid for general wavelengths. It will be shown 
that the effect of the anisotropic inertia is to damp the amplitude of the perturbation relative to the case where 
free-streaming neutrinos are absent, and that, in general, this damping depends on Q. Thus, it is important to 
have a solution capable of describing the damping for all wavelengths. For Q 2 3> 1, it has been shown that the 
solution for x( u ) can be written as a series of spherical Besscl functions whose coefficients are independent of 
Q [||. In [If, Weinberg analyzes the Q 2 3> 1 limit, provides results for the damping factor when Q 2 <C 1 and 
makes some observations about the damping forgeneral value of Q. Our aim here is to show that it is possible to 
extend the spherical Bessel function expansion p, 0] for x(u) to the case of general Q. Having done this, we can 
recapture the Q 2 < 1 results of [H , the Q 2 ^> 1 results of p and indicate how to obtain results for intermediate 
values of Q. 

In the next section, we derive the equation that must be satisfied by a spherical Bessel function expansion for 
examine its form in the Q 2 3> 1 and obtain a recurrence relation for general Q. In Section 3, we examine the 
long wavelength (Q 2 <g; 1) limit and determine the damping factor for this case. Section 4 contains a discussion 
of general wavelength case and we end with some conclusions. Lists of the various expansion coefficients and the 
details related to the solution of the recurrence relation are contained in the Appendix. 



2. SOLUTION OF THE GENERAL WAVELENGTH EQUATION 

In terms of the variable u, Eq. (J3j> becomes 

+ ^l±^91 x '(u) + X (u) = - T^ X2 f K(u u') X V) dv! , (6) 

with C = 24/„(0) = 9.72552 and the initial conditions 

X (0) = I, X'(0) = 0. (7) 

The kernel K is a linear combination of spherical Bessel functions 

1 f 1 j_ /P , 2\2 ixz sin(z) 3cos(2;) | 3sin(z) 1 . f ^ 2 . f ^ 1 



K{Z) = 16 U dx (1 ~ X ) CIXZ = " + = 15 J0(Z) + 21 J2(Z) + 35 J4(Z) ' (8) 

which suggests a solution to Eq. © exists in the form of a series of spherical Bessel functions. Indeed the work 
done in Ref. on the short wavelength limit shows that the convolution of a series of spherical Bessel functions 
with the kernel returns another series of spherical Besscl functions, so we look for a solution of the form 

oo 

X( U ) = ^2°<-njn(u). (9) 

71=0 

Multiplying Eq. ^ through by (u 2 + 4Qu) 2 , expanding all the terms, and dividing by u 4 yields 

8Q 16Q 2 \ /4 24Q 32Q 2 \ 

1 + — + — — ) X («) + ( - + — + ) X («) 

, 8Q 16Q 2 \ 16CQ 2 [ u , , , , 

+ H 1 5- XW = t— \ K{u- u)x {u)du . (10) 

The action of the differential operator on j n (u) gives: 

1 + 8Q + 16Q^\ £_ + /4 + 24Q + 32Q^\ ± + ( 1+ 8Q + 16Q^ 
u u 2 ) du 2 \u u 2 u 3 J du \ u u 2 

= -j [n(n + 1)16Q 2 + n(n + 2)8Qu + n(n + 3)u 2 ] j n (u) - — [2u 2 (u + 4Q)] j„+i(«) , (11) 



3n(u) 



3 



and Eq. (fTUj) becomes 

oo oo 

a» [n(n + 1)16Q 2 + n(n + 2)8Qu + n(n + 3)u 2 ] j n (u) - 2 ^ a n (u 3 + AQu 2 )j n+1 {u) 

n=Q 

= -16CQ 2 Y a n [ U K(u - u') du' . (12) 



n=0 



If we divide through by the highest power, the Bcsscl function recurrence relation (AS 10.2.18)[ 

jn{z) = jn-l(z) + jn+l(z) 

z ' (2n + 1) 

can be used to eliminate the powers of u. Dividing through by u 3 yields 



(13) 



^-n(n + 1) + ^n(n + 2) + -n(n + 3) 



n=0 ^ ' 



^ r dj n (u') J , 

——^a n K(u-u)——du. (14) 



2.1. Short Wavelength Limit 



Eq. (|13l) can now be used iteratively to eliminate the powers of u and obtain an expansion purely in spherical 
Bessel functions. However before we proceed, let us check what happens in the short wavelength limit (Q ^> 1). 
In this limit the terms with Q 2 dominate and Eq. (|14p becomes 



y~] ot n n{n + l)j n {u) = -C ^ a„ / K(u - u') n du' . 
n=0 n=0 *' U 

The derivative of j„(z) is given by (AS 10.2.19)01 

djn(z) nj„-i{z) - (n + l)j n+ i(^) 



and Eq. (fl5j) becomes 



oo oo „ 

E a » »(" + = -° E (2^i) y 

n=0 n=0 v ' " ,u 



(2n + 1) 



du'K(u - «') [nj n -i(u') - (n + l)j n+1 (u')] 



We can shift indices on the Bessel functions if we define a n — for n < to obtain 



oo oo / , u 

a n n(n + l)j n (u) = -C ^ u n J du' K(u- u')j n (u') , 



with 



a n+ i(n+l) noi n _i 



(2n + 3) (2n-l)_ 
The convolution integral can be evaluated using the technique given in Rcf. 



oo „ u oo 

^ loi I du' K(u- u') it (u') = ^ e « Jn ( u ) 5 

£=0 n=0 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 
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where 



2 £ 4^H)' + '" +1 C, , (2i) 



,£=0 m=0,2,4 



and 



with 



= y d»P„(s) [Qm(«)fll(*) + i , m(*)0<(*)] , 



do = i' d2 = ^' d4 = i- (23) 



2J a„ n(n + = -C 2J e « Jn(«) ■ (24) 

n=0 n=0 

This yields the following recurrence relation 

a n n(n+ 1) = -Ce n . (25) 



Then, Eq. ([18) becomes 



The e„ can be found using Eqs. (|2Tj). ([22]) and (|23|). The first 15 are given in Appendix A. For n = 0, Eq. (|25|) 
vanishes, so ao is undetermined. However, the initial condition x(0) = 1 fixes it to unity. The next six equations 
are 



2ai 



Cai 
"LIT 



6a 2 = -C f'g - ? 



12a 3 = -C* 

20a 4 = -C 

30a 5 = -C* 

42a 6 = -C 



15 6 

03 _ «i 

15 18 

10 20 15/ 

llai 33a3 as 
""225 700" + 15, 

13«o 13a2 143a4 a§\ 
~ 900 ~ "315 3150 + 15 ) 



The first equation requires ai = 0, and we see that all the odd terms depend recursively on ai. Thus, all the 
odd terms vanish. The constant C was defined as 

C = 24/„(0) = 9.72552. (26) 

Solving for the first 3 non-zero a n (n > 0) gives 

a 2 = 0.243807 

a 4 = 5.28424 x 10~ 2 

a 6 = 6.13545 x 10~ 3 



in agreement with Ref. . 
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2.2. Recurrence Relation 



Returning to the general case given by Eq. (fT4")l . we can apply Eq. (j2T))) to the righthand side and obtain 



n=0 



^-n(n + 1) + ^n(n + 2) + -n(n + 3) 



j n (u) - 2 y a n I 1 



n=0 



4Q 



2^ £ nJ«( M )- ( 27 ) 



n=0 



Eq. (fT3")) can then be used recursively to obtain (see Appendix C) 

oo oo 

[l6Q 2 /3„ + SQhn - n ) + S n - 2a n _ x ] j n (u) = -16CQ 2 K 3n{u) 



(28) 



where /3„, 7„, 6 n , 9 ni X n are defined in Appendix C and presented in Appendix D. Eq. (f2"5)l leads to the recurrence 
relation 



16Q 2 (/3„ + CX n ) + 8Q( 7 „ - 6 n ) + 8 n - 2a n -x = . 



(29) 



The equations given by Eq. (f2"9")) determine the coefficients of the spherical Bessel functions in Eq. © . With 
these coefficients, Eq. © provides an analytic solution to the inhomogeneous Eq. ^ for general values of the 
wave number Q. Eq. (|2"9"|) can also provide the coefficients of the homogeneous solution (the solution without 
free streaming neutrinos) by setting C = 0, yielding the following recurrence relation 



16Q 2 /?„ + 8Q( 7n - n ) + S n 



2a„_i = 0. 



(30) 



The coefficients were found by solving Eqs. (|29p and ppj) using Mathematica (with the condition a n = 
for n < 0). The first 15 inhomogeneous and homogeneous coefficients are given in Appendix A. The a n are the 
inhomogeneous coefficients and the b n are the homogeneous coefficients. This convention will be used throughout 
the remainder of the paper. 



3. THE LONG WAVELENGTH LIMIT 



We found an analytic solution to Eq. © consisting of an expansion in spherical Bessel functions of the following 
form 

oo 

X(u) = ^2a n (Q)j n (u) , (31) 

where the a n (Q) are explicit functions of the wave number Q. The convention used here will be to replace a n (Q) 
with a n (Q) for the inhomogeneous solution and b n {Q) for the homogeneous solution. Upon first glance, this 
expansion seems to be divergent in the limit Q — > since the coefficients depend inversely on Q. However, the 
expansion is in fact finite due to the implicit presence of Q in j n (u). To make the low Q limit unambiguous, the 
Q dependence will be removed from the Bessel functions using the following theorem (AS 9.1.74)Q 

\-"M\z) =J2 TO , (2) W*)> ( 32 ) 

m=0 

where fi and A can be arbitrary complex numbers. If we define u = Qs, where s = 2 (^Jl + y — l) , then we can 
write Eq. (f3"Tj) as 



00 

x(s,Q) = ^2a n (Q)j„(Qs) 

71=0 



(33) 
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Using Eq. (|32p we can write 



(1-Q 2 ) m ( s 



(|) 3n+m{s) 



(34) 



so Eq. ([55)1 becomes 



^,Q) = X>„(Q)Q"]T 



(i-Q 2 ) 



2\rn 



n=0 



rn—0 



(|) Jn+m(s) • 



(35) 



It is now convenient to define new coefficients c n (Q) (inhomogencous) and d n (Q) (homogeneous) that arc finite 
as Q -> 



c n (Q) = a n {Q)Q n 
d n {Q) = b n (Q)Q n 



(36) 
(37) 



The first 15 c n and d n are given in Appendix B. The expansion in Eq. (|35[) becomes (in the inhomogencous case) 

(i-Q 2 r 



c(*,g) = X; c »(o)E 



n=0 



m=0 



2 J in+m(s) 



(38) 



For Q — ;> 0, only c contributes and the canonical damping factor \x( s l)/xo( s l)\ 2 —> 1 since the homogeneous 
expansion is the same in this limit (i.e. only do contributes and it is also 1). We would also like to examine 
the Weinberg damping factor \x' '(s l) / Xo( s l)\ 2 and show that it is well defined and finite in the limit Q — > 0. 
Differentiating Eq. (|38[) with respect to s yields 



x'( S ,Q) = ^c„(Q)^ 



n=0 



m=0 



(i-Q 2 ) T 

2 m m! 



ms m - i j„ +m (s)+s^ 



(n + m)j ra+m _i(3) - (n + m + l)j n+m+ i(s) 
2(n + m) + 1 



(39) 



The cosmologically interesting value of s is when y = 22.1 fljf/i 2 . This is the value of s at last scattering, denoted 
Sl,. For il]\jh 2 = 0.15 we have = 2.15452. Plugging this value of s into Eq. ([3^]) and expanding using 50 terms 
in each sum yields 

x(s Ll Q) = -0.573659 Q 2 + 0.243294 Q 4 - 0.0381658 Q 6 + 0.00315913 Q 8 + 0{Q W ) . (40) 
For Q -C 1, the Q 2 term dominates and we have 

x'(s L ,Q) pa -0.573659 Q 2 . (41) 

The c n (Q) in Eq. Q39[) can be replaced with the d n (Q) to obtain an expansion of the same type for the 
homogeneous solution. Using the same method, we obtain the following expansion 

Xo(sl, Q) = -0.601268 Q 2 + 0.264497 Q 4 - 0.0424247 Q 6 + 0.00356526 Q 8 + 0(Q 10 ) . (42) 

As before, for Q -c 1, the Q 2 term dominates and we have 

Xo(sl 7 Q) « -0.601268 Q 2 (43) 

so the Weinberg damping factor |x'( s i; Q)/x'o( s l, Q)\ 2 for Q = can be computed as 



x'(s L ,0) 



lim 



-0.573659 Q s 



-0.601268 Q 2 



0.910273. 



(44) 
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in agreement with Weinberg's discussion on the low Q limit. As Weinberg remarks "This damping is relatively 
insensitive to Q for small Q." Confirmation of this statement is seen here where the damping is independent of 
Q until the fourth order terms become important. The finiteness of this damping factor shows that the solution 
is well defined for arbitrarily small values of Q. 

The expansions given in Eqs. (|38|) and (f39|) are valid for general values of Q. However for Q > 1 the expansion 
quickly starts to fail with any finite number of terms, due to the fact that the large powers of Q overpower their 
coefficients and the final term in any finite sum dominates. Writing out Eqs. (|38|) and (|39[) using Mathematica 
(100 terms per sum) will give an expansion in powers of Q up to Q 300 . As shown in Figurc[TJ this expansion is 
well behaved for Q < 1 and converges up to approximately Q = 15. 




FIG. 1: xoO-L,Q)(blue) and X (sl, Q)(red). 

The ratios of x' s an d their derivatives can be used to plot the damping factors, as shown in Figure[2] These 
damping factors are the ratios of two oscillating functions and therefore become infinite at the zeros of the function 
in the denominator. However, since the functions arc only slightly out of phase, these spikes occur where x'{ s l) 
is small and according to Ref. this makes the spikes uninteresting since the multipole coefficients for the 
corresponding values of £ will be very difficult to measure. 





FIG. 2: The damping factor \x(sl)/xo{sl)\ 2 is shown in the left panel and \x' ( s l)/x'o{ s l)\ 2 is shown in the right panel. 
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\X\ s L)/X0(SL)\ 


IX ( s l)/Xo( s l )\ 


n 

u 


1 

1 


u.yiuz i o 


0.01 


1.00001 


0.910275 


0.1 


1.00081 


0.910561 


0.55 


1.02745 


0.919509 


0.8 


1.06802 


0.931291 


1 


1.13058 


0.946036 


5 


0.669293 


1.23524 


10 


0.811905 


0.644909 



TABLE I: Damping values for various Q's from Eq. (|38[) are listed. 



4. DAMPING FOR GENERAL WAVELENGTHS 



Though the expansion given in Eq. (|38p is valid for general values of Q, it is cumbersome to use in practice. 
However, now that the solution has been shown to be finite in the long wavelength limit, we can return to a 
more useful expansion of the form given in Eq. (|33l) , 

oo 

X (.s,Q) = 5>„(Q)j„(Q S ). (45) 

?i=0 

Here s = 2 (-y/l + y — l) and the a n (Q) should be replaced with the a n (Q) (b n (Q)) for the inhomogeneous 
(homogeneous) solution. As long as numerical applications of Eq. (|45|) explicitly include the Q in the argument 
of the Bessel functions, this expansion reproduces the damping found in the previous section (given in TableU]) 
as Q — > 0, despite the fact that the coefficients diverge. This is due to the fact that j n (Qs) goes like Q n for small 
Q, removing the divergences produced by the coefficients. For Q > 1, there is no difficulty and asymptotically 
this expansion reproduces the damping found in Ref. @ for Q 3> 1, as shown in Figure[3l 

Unless otherwise specified, 100 terms were used in all numerical evaluations of Eq. (|45[) occurring in this paper. 
The observation in Ref. [IJ that the reduction in \x' ( s l) / Xo( s l)\ 2 from 1 for Q = 0.55 and Q = 0.8 is about 
8% and 7% respectively is confirmed in TablclLTl Also shown in this Table is the trend of \x' ( s l)/x'o( s l)\ 2 in the 
relatively flat regions between the spikes steadily decrease from the value « 0.9 for Q <C 1 to a value close to 
0.644 for Q w 10. 
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FIG. 3: The damping factor |x'( s i)/Xo( s -t)| 2 is shown along with the horizontal line representing A 2 = 0.645019. 



9 





\X( s L)/X0(SL)\ 




U.Ul 


i nnnni 

l.UUUUl 


n Qi no7^ 


U. 1 


l.UUUol 


u.y iuooi 


u.oo 


l.UZ / 40 




u.o 


i nfisn? 

l.uUOUi 


n Q^l 9Q1 


1 


1.13058 


0.946036 


10 


0.811905 


0.644909 


10 2 


0.64902 


0.866294 


10 3 


0.646315 


0.642628 


10 4 


0.644227 


0.64506 


10 5 


0.645004 


0.645055 


10 6 


0.645013 


0.645008 



TABLE II: Damping values for various Q's from Eq. (|45|l 
5. CONCLUSIONS 

We have shown that the treatment of gravitational wave damping by free streaming neutrinos can be framed 
in terms of a series of spherical Besscl functions for all values of the reduced wave number Q. The result for the 
coefficients of the spherical Bessel series when Q ^> 1 Q emerges quite simply from the coefficient recurrence 
relation for a general Q. For the opposite limit, Q <C 1, an identity is used to remove the <5™-dependence from 
jn(Qs) and the analysis can then be extended to arbitrarily small values of Q, including the non-leading behavior 
of the series. Once Q is large enough, the coefficients obtained from the general recurrence relation can be used 
directly, even though they contain inverse powers of Q. 
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Appendix A: Expansion Coefficients 







The relations between the e n (Q) and the expansion coefficients a n (Q) are listed below for n = to 14. 
eo(Q) 
ei(Q) 

ea(Q) 
e 3 (Q) 
e 4 (<9) 
e 5 (<9) 
ee(<9) 
er(Q) 
es(Q) 
eg(Q) 
eio(Q) 

cii(Q) 

ei 2 (<9) 

cis(g) 
£i 4 (g) 









15 






a 2 (Q) 


ao(<9) 




15 


6 






a x {Q) 




15 


18 




a 4 (Q) 


a,2(Q) Oo(Q) 


15 




1 n 




33a 3 (Q) 


lloi(Q) 


15 


700 


225 




143a 4 (Q) 


13a 2 (Q) 


15 


3150 


315 


a 7 (Q) 


13a 5 (Q) 


lla 3 (Q) 


15 


294 


294 


an(Q) 


17a 6 (Q) 


221a 4 (Q) 


15 


392 


6300 


aa(Q) 


323a 7 (Q) 


19a 5 (Q) 


15 


7560 


567 


aio(Q) 


19a 8 (Q) 


17a 6 (Q) 


15 


450 


525 


an(Q) 


23a 9 (Q) 


437a 7 (Q) 


15 


550 


13860 


ai2(Q) 


ii5o 10 (g) 


5a 8 (Q) 


15 


2772 


162 


ai3(Q) 


15ou(Q) 


69a 9 (Q) 


15 


364 


2275 


au{Q) 


261a 12 (Q) 


145oi ( 


15 


6370 


4851 



13a (Q) 
900 
5ai(Q) 
588 

17o 2 (Q) 
2520 T " 

4693a 3 (g) 
793800 

41a 4 (Q) 
7560 



17ao(Q) 
11025 
19oi(Q) 
T 18900 
a 2 (Q) a (Q) 



1260 



2520 



69989a 5 (Q) 299a 3 (Q) 23a x (Q) 



13721400 
18715a 6 (g) 

3841992 
887a 7 (Q) 



436590 
85a 4 (Q) 
^ 137214 " 
13319a 5 (Q) 



83160 
a 2 (Q) 25a (Q) 
4536 174636 
1751a 3 (Q) ai(Q) 



188760 23123100 9249240 9555 
2850091a 8 (Q) 8207a 6 (g) 551a 4 (Q) 



624323700 



15057900 



3243240 



+ 



29a 2 (Q) 29a (0) 



343035 



463320 
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Using C = 24/„(0) = 9.72552, the first fifteen a n (Q) are 

ao(Q) = 1 
oi (Q) = 
a 2 (Q) = 0.243807 



a 3 (Q) 



4 (<9) = 0.0528424 

05 (<9) 



Og(<9) 



on(Q) 



013(0) 



0.843856 
Q 

1.31658 



Q 2 

2.60385 1.98354 



a 6 (Q) = 0.00613545 
or(Q) 



Q 3 

6.16442 8.15545 



Q 4 Q 2 
16.7669 31.6633 2.81654 



Q 5 Q 3 Q 

50.564 130.451 22.4418 

a 8 (Q) = 0.000297534- 



Q 6 Q 4 Q 2 

162.551 583.689 141.989 3.66375 



Q 7 Q 5 Q 3 Q 

524.266 2845.09 868.128 46.8401 

aio(<9) = 0.0000616273- 



o 8 o 6 o 4 o 2 

1458.97 15066.7 5429.37 433.843 4.49219 



Q 9 Q 7 Q 5 Q 3 Q 



, ™™™ 1012.45 86302.4 35476. 3653.68 83.7481 

ai 2 (<9) = -4.998662 x 10" 6 - 



Q 10 Q s Q 6 Q 4 
39324.1 532143. 244196. 30143.3 1064.03 5.3138 



a 14 (Q) = 2.33661 x 10 



11 Q 9 Q 7 Q 5 Q 3 Q 

570180. 3.5156 x 10 6 1.7758 x 10 6 251558. 



-6 



Q 12 Q 10 Q 8 Q 6 

11768.6 135.604 



Q 4 Q 2 

For large Q, the coefficients of the a 2n ((3) reduce to the results in Ref. |f|. For a fixed value of y, these coefficients 
and the b n (Q) below yield finite results. This is illustrated in Figure |H 
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The coefficients for the homogeneous solution are 



HQ) 

h(Q) 
b2(Q) 

HQ) 
64(g) 
h(Q) 
b 6 (Q) 
67(g) 

HQ) 
HQ) 

610(G) 
b n (Q) 
612(g) 
613(g) 
614(g) 



1 



= 



1.45833 

~~ Q~ 

2.95313 
7.57969 
23.4609 



84.8364 

~Q^~ 

350.539 



2.56667 

~~ Q~ 

12.4583 

+ _ w~ 

55.7289 

g 3 + 

261.588 



1628.06 1323.89 



Q 7 

8391.91 

47522.4 

~Q^~ 
293207. 

g^ + " 

1.9574 x 10 6 



1.4056 x 10 7 



g 5 

7261.36 



g 6 

43105.3 



Q 7 

275999. 



3.61607 

~~ Q~ 

31.5725 
218.214 



1451.12 
9831.73 



69352.2 



Q s 



Q 6 



1.89854 x 10 6 



+ 



4.64444 
63.3263 

+ — 

624.673 

— + 

5591.11 

+ — 

513799. 



5.66288 



110.739 



48915.7 1475.69 6.67582 



1.39746 x 10 7 



Q 



10 



g 7 g 5 

4.01079 x 10 6 



g 3 

431978. 
~Q^' 



17104.2 176.825 



g 2 




FIG. 4: The plots for x( s l,Q) (dashed) and Xo(sl,Q) (solid) are shown using the 15 coefficients listed above. 
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Appendix B: Long Wavelength Finite Expansion Coefficients 



The coefficients for the inhomogeneous solution are 



c (Q) 


= 1 




ci(Q) 


= 




„ /n\ 

C2{Q) 


= 0.243807Q 2 




= 0.843856Q 2 


cAQ) 


= Q 2 


(0.0528424Q 2 - 1.31658) 




= Q 2 


(2.60385- 1.98354Q 2 ) 


MQ) 


= Q 2 


(0.00613545Q 4 + 8.15545Q 2 - 6.16442) 


cr(Q) 


= Q 2 


(2.81654Q 4 - 31.6633Q 2 + 16.7669) 


MQ) 


= Q 2 


(0.000297534Q 6 - 22.4418Q 4 + 130.451Q 2 - 50.564) 


cb(Q) 


= Q 2 


(-3.66375Q 6 + 141.989Q 4 - 583.689Q 2 + 162.551) 


cio(Q) 


= Q 2 


(0.0000616273Q 8 + 46.8401Q 6 - 868.128Q 4 + 2845.09Q 2 - 524.266) 


cn(Q) 


= Q 2 


(4.49219Q 8 - 433.843Q 6 + 5429.37Q 4 - 15066. 7Q 2 + 1458.97) 


cu(Q) 


= Q 2 


(-4.99866 x 10- 6 Q 10 - 83.7481Q 8 + 3653.68Q 6 - 35476Q 4 + 86302.4Q 2 - 


cis(Q) 


= Q 2 


(-5.3138Q 10 + 1064.03Q 8 - 30143.3Q 6 + 244196Q 4 - 532143Q 2 - 39324.1 


Cl4(<3) 


= Q 2 


(2.33661 x 10~ 6 Q 12 + 135.604Q 10 - 11768.6Q 8 + 251558Q 6 - 1.7758 x 10 6 



1012.45) 



570180) 



The coefficients for the homogeneous solution are 



MQ) 


= 1 






di(Q) 


= 






d 2 {Q) 


= 






MQ) 


= 1.45833Q 2 






MQ) 


= -2.95313Q 2 






MQ) 


= Q 2 (7.57969- 2.56667Q 2 ) 






MQ) 


= Q 2 (12.4583Q 2 - 23.4609) 






MQ) 


= Q 2 (3.61607Q 4 - 55.7289Q 2 + 


84.8364) 




MQ) 


= Q 2 (-31.5725Q 4 + 261.588Q 2 


- 350.539) 




MQ) 


= Q 2 (-4.64444Q 6 + 218.214Q 4 


- 1323.89Q 2 


+ 1628.06) 


d w (Q) 


= Q 2 (63.3263Q 6 - 1451.12Q 4 + 


7261.36Q 2 - 


8391.91) 


du(Q) 


= Q 2 (5.66288Q 8 - 624.673Q 6 + 


9831. 73Q 4 - 


43105.3Q 2 + 47522.4) 


di 2 {Q) 


= Q 2 (-110.739Q 8 + 5591. HQ 6 


- 69352.2Q 4 


+ 275999Q 2 - 293207) 


di 3 (Q) 


= Q 2 (-6.67582g 10 + 1475.690* 


! - 48915. 7Q ( 


+ 513799Q 4 - 1.89854 x 10 6 Q 2 + 1.9574 x 10 6 ) 


d u {Q) 


= Q 2 (I76.825g 10 - 17104.2Q 8 - 


h 431978Q 6 - 


4.01079 x 10 6 g 4 + 1.39746 x 10 7 Q 2 - 1.4056 x 10 7 



In contrast to Figure [T] using this set of 15 coefficients in Eq. (|38|) gives a convergent expansion to about Q = 5, 
as illustrated in Figure [5] This is adequate to treat the Q < 1 region. 
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FIG. 5: The plots for x{ s l, Q) (dashed) and Xo(sl,Q) (solid) are shown using the 15 coefficients listed above in Eq. (|38[) . 



Appendix C: Recurrence Calculation 



We want to eliminate powers of u from 

oo 



n=0 



l^! n (n + 1) + ^-n(n + 2) + -n(n + 3) 

vr u u 



n=0 



( u )-2*£a n (l + ^-)j n+1 ( 



16CQ 2 



X>„j„(u) (CI) 



n=Q 



The Bcssel function recurrence relation 

jn(z) _ jn-l(z) + jn+l(z) 

z ~ (2n+l) 
can be used recursively to obtain the following relations: 

3n(z) 3n-l{z) 2jn(z) 



3n+2(z) 



z 2 (2n- l)(2n + 1) (2n - l)(2n + 3) (2n + l)(2n + 3) 

jn(f) _ jn-3(^) 3jn-l{z) 

z 3 ~ (2n + l)(2n - l)(2n -3) (2n - 3)(2n + l)(2n + 3) 

3jn+i(z) j»+3(z) 

(2n- l)(2n + l)(2n + 5) (2n + 1) (2n + 3) (2n + 5) 

Consider the first term on the LHS of Eq. (|C1[) 



(C2) 



(C3) 



(C4) 



16Q 2 , N sg , , 1 . 

— f-n(n + 1 + -|-n(n + 2) + -n(n + 3 



in(w) 



16Q 2 ^ o„ n(n + 1) + 8Q £ a n n{n + 2) + ]T a„ n(n + 3) 



3n(u) 



3n{u) 



(C5) 



n— n— n— 

The first term in Eq. (|C5|) can be re- written as follows with the aid of Eq. (|C4[) : 



16Q 2 ]T a„ n(n + 1) = 16Q 2 a n n(n + 1) 



n=0 



n=0 



3n-3(z) 



3jn-l(z) 



(2n + l)(2n - l)(2n - 3) (2n - 3)(2n + l)(2n + 3) 

3jn+l(z) in+3(z) 



(2n - l)(2n + l)(2n + 5) (2n + l)(2n + 3)(2n + 5) 



(C6) 
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If we define a n = for n < then the indices on the Bessel functions can be shifted to obtain: 



16Q 2 «n n(n + 1) ^ = 16Q 2 £ j„ («) 



n=0 



So we can write: 



n=0 



a n+3 (n + 3)(n + 4) 3 a n+i (n + l)(n + 2) 



+ 



(2n + 3)(2n + 5)(2n + 7) (2n - l)(2n + 3)(2n + 5) 
3a n -in(n — 1) a„_ 3 (n - 3)(n - 2) 



(2n-3)(2n- l)(2n + 3) (2n - 5)(2n - 3)(2n - 1) 



16Q 2 £ a„ n(n + 1) = 16 ^ 2 Pnjn{u) 



n=0 



n=0 



with: 



a„ +3 (n + 3)(n + 4) 3 a„+i (n + l)(n + 2) 



(2n + 3)(2n + 5)(2n + 7) (2n - l)(2n + 3)(2n + 5) 
3a„_in(n — 1) a„_3 (n - 3) (re - 2) 



(2n-3)(2n-l)(2n + 3) (2n - 5)(2n - 3)(2n - 1)_ 
Following the same procedure for the rest of Eq. (|C5|) using Eqs. (|C2[) and (|C3|) yields: 



n=0 



16Q 5 



1 



n(n + 1) + —kn(n + 2) + -n(n + 3) 



n=0 



with: 



a n+2 (n + 2)(n + 4) 2a„n(n + 2) a„_ 2 n(n — 2) 



(2?i + 3)(2n + 5) (2n- l)(2n + 3) (2n-3)(2n-l) 

, _ [ ctn+i (n + l)(n + 4) a ra _i (n- l)(n + 2) 
n_ [ (2n + 3) (2n-l) 

Similarly, the second term on the LHS of Eq. (|C1|) becomes: 

4Q X 



n=0 



2 a n 1 + — Jn+1 («) = -2 ^ [4Q0„ 



+ Q!n-1 



with: 



n=0 



«n-2 



(2n + 3) (2n-l) 



So Eq. (|C1|) becomes: 



^ [l6Q 2 /3„ + 8Q( 7n - On) + $n ~ 2a„_i] j n (u) 



16CQ 2 



n=0 



n=0 



(C7) 



(C8) 



(C9) 



= [i6Q 2 /?„ + 8g 7 „ + <y n ] (cio) 



(Cll) 
(C12) 

(C13) 
(C14) 

(C15) 



With the condition a n = for n < 0. The e„ are already zero for n < 0, so the same method as above can be 
applied to the RHS of Eq. (|CT5]) . Using Eq. {CD gives: 



16CQ 2 



]T}e„j n (u) = -16CQ 2 j„(?l) 



n=0 



n=0 



£n+3 



3 e n +i 



_ (2n + 3)(2n + 5)(2n + 7) (2n - l)(2n + 3)(2n + 5) 
3e„_i e n -3 



(2n - 3)(2n - l)(2n + 3) (2n - 5)(2n - 3)(2n - 1) 



(C16) 



So the RHS of Eq. (ICTBT) becomes: 

16CQ 2 



^ enjn(u) = -16CQ 2 \ n jn{u) 



with: 



n=0 



e n+3 



n=0 



3e n 



(2n + 3)(2n + 5)(2n + 7) (2n - l)(2n + 3)(2n + 5) 
3e„_i e n -3 



(2n - 3)(2n - l)(2n + 3) (2n - 5)(2n - 3)(2n - 1)_ 
Substituting, Eq. (|C15|) becomes 

oo oo 

[l6Q 2 /?„ + 8Q( 7 „ - 9 n ) + S n - 2a n _i] j n (u) = -16CQ 2 ^ A„ j„( 



n=0 



71 = 



which is Eq. (|28|) of this paper. 



Appendix D: Recurrence Coefficients 



a„ +3 (n + 3)(n + 4) 3 a n+1 (n + l)(n + 2) 



(2n + 3)(2n + 5)(2n + 7) (2n - l)(2n + 3)(2n + 5) 
3a„_i n(n - 1) a„_ 3 (n - 3)(n - 2) 



(2n - 3)(2n - l)(2ra + 3) (2n - 5)(2n - 3)(2n - 1) 



7n 



a«+2 (n + 2)(n + 4) 2a n n(n + 2) a n -2 n.(n — 2) 



(2n + 3)(2n + 5) (2n-l)(2n + 3) (2n - 3) (2n - 1) 



a„ + i (n+ l)(n + 4) a n _i (n - l)(n + 2) 



(2n + 3) 



(277.- 1) 



a„ a„_2 



(2n + 3) (2n-l) 



£n+3 



3e„ 



{2n + 3)(2?i + 5)(2n + 7) (2n - l)(2n + 3)(2n + 5) 
3e„_i e n -3 



(2n - 3)(2n - l)(2n + 3) (2n - 5){2n - 3)(2n - 1) 



